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Abstract 

The  deformation  of  rotorcraft  blades  are  frequently  clas¬ 
sified  as  moderate,  exceeding  the  small  deformation  lim¬ 
itation  of  typical  beam  finite  elements.  Advances  in 
computer  capabilities  have  made  possible  an  implemen¬ 
tation  of  a  hybrid  element  combining  rigid  and  flexible 
body  kinematics  in  the  Rotorcraft  Comprehensive  Anal¬ 
ysis  System  (RCAS).  Two  challenges  to  new  software 
are  to  verify  the  correctness  of  the  implementation  and 
to  ascertain  whether  it  represents  an  improvement.  Com¬ 
parisons  are  made  to  analytical  results,  to  measurements 
from  the  “Princeton  Beam”  and  Maryland  vacuum  cham¬ 
ber  experiments,  and  to  UH-60  flight  test  data.  Inspection 
of  these  comparisons  lead  to  the  conclusion  that  the  re¬ 
sults  obtained  with  RCAS  correlate  well  with  analytical 
and  experimental  results,  sometimes  representing  signif¬ 
icant  improvements. 

Introduction 

The  broad  discipline  of  rotorcraft  dynamics  involves  a 
variety  of  important  time  dependent  and  nonlinear  phe¬ 
nomena.  Even  when  the  major  challenges  of  aerody¬ 
namics  and  material  characteristics  are  ignored,  there  are 
many  complex  fundamental  problems  involving  purely 
mechanics  and  structural  dynamics  phenomena,  includ¬ 
ing  rotating  systems,  large-motion  multi-body  dynam¬ 
ics,  and  moderate-to-large  deformations  of  elastic  beams 
of  arbitrary  shape.  In  the  past  these  types  of  problems 
have  often  constituted  a  significant  challenge  to  success¬ 
ful  and  accurate  treatment  of  important  rotorcraft  dynam¬ 
ics  and  aeroelastic  phenomena.  With  the  continuing  de- 
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velopment  of  advanced  comprehensive  analysis  methods, 
many  of  these  problems  are  now  amenable  to  rigorous 
mathematical  treatment  and  this  provides  new  capabili¬ 
ties  for  treating  many  rotorcraft  dynamics  problems  — 
although  not  all  of  the  challenges  have  been  reduced  to 
the  category  of  routine  engineering. 

This  paper  will  examine  a  number  of  selected  problems 
relevant  to  analysis  of  rotorcraft  structural  mechanics  and 
dynamics  that  deal  with  the  issues  mentioned  above.  The 
focus  of  this  examination  will  encompass  three  related 
themes: 

1.  a  consideration  of  the  methods  and  philosophy 
of  demonstrating  and  validating  rotorcraft  analysis 
computer  programs,  and 

2.  a  careful  re-examination  of  a  number  of  important 
fundamental  laboratory  experiments, 

3.  a  presentation  of  results  providing  initial  validation 
of  the  RCAS  computer  program. 

This  introductory  section  first  provides  the  background 
leading  to  the  development  of  a  computer  program  which 
implements  a  hybrid  finite  element  encompassing  both 
rigid  and  flexible  body  kinematics.  It  then  describes  the 
implementation  of  the  hybrid  element  in  a  computer  pro¬ 
gram.  Next,  it  identifies  the  methods  used  for  demonstrat¬ 
ing  and  validating  the  element  formulation  and  software 
implementation.  Finally,  it  describes  the  representative 
problems  which  were  selected  for  this  paper. 

This  introductory  section  is  followed  by  four  sections 
of  results  comparisons.  The  sections  treat  groups  of  prob¬ 
lems  of  the  following  four  types:  static  elasticity  prob¬ 
lems,  fixed  reference  frame  eigenvalue  problems,  rotat¬ 
ing  reference  frame  eigenvalue  problems,  and  periodic 
boundary  load  problems. 

Conclusions  follow  the  four  results  comparison  sec¬ 
tions. 
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Background 

The  presence  of  articulated  flexible  structures  such  as 
solar  arrays,  booms,  and  antennae  on  spacecraft  stim¬ 
ulated  the  development  of  computational  methods  cou¬ 
pling  rigid  and  flexible  body  dynamics.  Early  methods  al¬ 
lowed  articulated  central  rigid  bodies  (typically  restricted 
to  a  tree  topology)  with  flexible  appendages.  Ultimately, 
finite  element  and  modal  analysis  methods  were  permit¬ 
ted  for  the  appendages  [1]. 

The  governing  equations  were  subsequently  reformu¬ 
lated  to  incorporate  rigid  body  and  flexible  body  degrees 
of  freedom  in  each  body,  and  to  allow  arbitrary  intercon¬ 
nection  topologies  [2],  This  methodology  was  incorpo¬ 
rated  in  the  General  Rotorcraft  Aeromechanical  Stability 
Program  (GRASP)  [3,  4],  GRASP  computed  eigensolu- 
tions  for  a  rotorcraft  in  vertical  flight  or  ground  contact. 
The  run  times  required  to  calculate  the  motion  of  the  ref¬ 
erence  frame  degrees  of  freedom  and  the  flexible  body 
degrees  of  freedom  associated  with  each  element  were 
high,  even  on  the  Cray  supercomputer  for  which  it  was 
written. 

Using  earlier  technology,  the  comprehensive  Second 
Generation  Comprehensive  Helicopter  Analysis  System 
(2GCHAS)  [5,  6,  7,  8]  ran  faster  than  GRASP,  but  run 
times  remained  a  concern.  However,  it  had  difficulties 
with  some  rigid  body  motions  and  moderate  deforma¬ 
tions.  The  displacement  fields  which  can  be  represented 
by  2GCHAS  (and  other  programs)  inadequately  repre¬ 
sented  the  combination  of  rigid  body  motions  of  assem¬ 
blages  of  elements  experiencing  moderate  deformations. 

The  Rotorcraft  Comprehensive  Analysis  System 
(RCAS)  [9,  10,  11]  is  an  Army  sponsored  computer 
program  developed  for  interdisciplinary  engineering,  re¬ 
search,  and  design  of  arbitrary  rotorcraft  configurations. 
Applications  include  performance,  dynamics,  aerody¬ 
namics,  aeroelastic  stability,  loads  and  vibration,  sta¬ 
bility  and  control,  and  acoustic  characteristics  of  rotor¬ 
craft.  RCAS  was  derived  from  the  Army’s  previously 
developed  Second  Generation  Comprehensive  Helicopter 
Analysis  System  (2GCHAS)  and  involved  significantly 
changing  important  parts  of  the  mathematical  basis  and 
the  system  architecture  of  the  2GCHAS  system.  These 
revisions  were  designed  to  overcome  specific  technical 
limitations  of  the  predecessor  2GCHAS  system  and  have 
resulted  in  substantial  improvements  to  the  computa¬ 
tional  efficiency  of  the  system  in  addition  to  significantly 
expanding  the  functionality  and  robustness  of  the  system. 


Implementation 

The  principal  changes  represented  in  the  RCAS  mathe¬ 
matical  basis  and  system  architecture  include  a  reformu¬ 
lation  of  the  elements  that  make  up  the  RCAS  element 
library  and  a  revised  scheme  in  which  finite  elements  are 
assembled  into  the  rotorcraft  structure,  as  described  in 
reference  [12].  The  new  element  associates  a  reference 
frame  a  nonlinear  beam  element  with  the  root  of  the  beam 
element  attached  to  the  origin  of  the  reference  frame.  Im¬ 
provements  in  computer  hardware  allow  rapid  calculation 
of  reference  frame  motion  as  well  as  flexible  body  defor¬ 
mations  for  every  element. 

Revised  algorithms  were  required  for  element  assem¬ 
bly.  Rather  than  assembling  together  all  blade  nonlinear 
beam  elements  referred  to  a  single  rotor  blade  reference 
frame,  the  revised  approach  allows  large  deformations  of 
the  blade  while  each  beam  element  undergoes  only  mod¬ 
erate  elastic  deformations  with  respect  to  its  own  refer¬ 
ence  frame.  Each  reference  frame  is  attached  to  the  tip 
extremity  of  the  adjacent  beam  element. 

An  important  benefit  of  the  additional  reference  frames 
and  the  revised  assembly  process  is  the  accuracy  and  ro¬ 
bustness  of  analyses  involving  large  displacements  of  the 
vehicle  rigid-body  degrees  of  freedom.  In  contrast,  only  a 
limited  capability  for  large  motion  maneuvers  is  available 
with  the  basic  2GCHAS  formulation.  It  also  makes  it  fea¬ 
sible  to  implement  a  practical  modal  reduction  approach 
for  nonlinear  rotor  blades  that  enables  significant  savings 
in  computational  time  for  typical  rotorcraft  problems. 

Validating  Finite  Elements 

There  are  two  fundamental  questions  that  must  be  an¬ 
swered  when  software  for  a  new  finite  element  is  devel¬ 
oped: 

1.  Does  the  software  correctly  implement  the  theory 
for  the  new  element,  the  assembly  of  elements,  and 
solution  of  the  resulting  equations? 

2.  Does  the  new  element  demonstrate  better  results 
than  previous  elements? 

The  theory  for  a  finite  element  is  usually  derived  based 
on  an  energy  formulation.  The  types  of  energy  terms  in¬ 
volved  are  elastic  energy,  kinetic  energy  (associated  both 
with  motion  relative  to  the  reference  frame  and  with  ref¬ 
erence  frame  motion),  and  the  work  done  (associated 
both  with  body  loads  and  boundary  loads). 

One  way  to  organize  the  addressing  of  the  above  ques¬ 
tions  is  to  answer  them  for  each  of  the  energy  terms  at 
the  element  level.  The  same  questions  also  need  to  be  an¬ 
swered  for  the  assembly  of  elements  discretizing  a  con¬ 
tinuum  (e.g.,  a  blade)  and  for  the  heterogeneous  assem- 


blages  of  various  elements  (e.g.,  a  lumped  mass)  and  sub- 
assemblies  (e.g.,  a  fuselage)  which  comprise  a  model. 

Some  of  these  questions  can  be  answered  through  ap¬ 
propriate  module  testing.  But,  ultimately,  the  demonstra¬ 
tion  and  validation  must  be  accomplished  by  comparison 
of  computed  results  with  known  results.  There  are  three 
types  of  sources  for  results: 

1 .  Analytical  results  are  available  for  some  simple  con¬ 
figurations.  Results  may  be  in  closed  form,  in  series 
form,  or  in  transcendental  form  with  reported  solu¬ 
tions.  Program  results  for  the  configurations  it  is 
designed  to  address  should  show  close  agreement. 

2.  Laboratory  experimental  results  obtained  under 
carefully  controlled  conditions  are  available  for 
more  complex,  but  still  relatively  basic,  configura¬ 
tions.  Program  results  should  fall  within  experi¬ 
mental  error  if  systematic  apparatus  errors  are  elim¬ 
inated. 

3.  Wind  tunnel  or  flight  test  results  are  available  for 
real  configurations.  The  complexity  of  real  sys¬ 
tems,  limited  instrumentation,  experimental  errors 
due  to  uncontrolled  test  conditions,  unmodeled  fea¬ 
tures,  and  uncertainties  in  system  properties  are  all 
likely  to  give  rise  to  discrepancies  with  program  re¬ 
sults.  However,  improvements  in  correlation  here 
are  the  acid  test  of  the  value  of  new  elements  and 
software. 

Selected  Problems 

Results  from  a  variety  of  analytical,  laboratory  and  field 
systems  were  selected  to  demonstrate  and  validate  the 
new  element  and  the  RCAS  program  which  implements 
it.  In  each  case,  RCAS  was  used  to  solve  the  compu¬ 
tational  problem  of  representing  the  system  and  predict¬ 
ing  the  resulting  behavior.  RCAS  computational  results 
were  then  compared  with  the  known  results  for  the  sys¬ 
tem.  One  or  more  input  files  specifying  the  RCAS  com¬ 
putational  model  were  written,  RCAS  was  run,  and  the 
results  were  extracted  from  graphics  related  files.  The 
representative  comparisons  presented  here  are  grouped 
by  the  finite  element  energy  terms  which  the  problems 
exercise. 

Only  elastic  energy  appears  in  the  first  group  of  prob¬ 
lems.  Their  solutions  are  static  deformations.  Compar¬ 
isons  are  made  to  static  “Princeton  Beam”  experiment 
data  and  two  analytical  elastica  problems. 

The  second  group  of  problems  adds  kinetic  energy 
terms  for  motion  relative  to  a  fixed  reference  frame.  The 
solutions  are  natural  frequencies.  Comparisons  are  made 
to  dynamic  “Princeton  Beam”  experiment  data  and  an  an¬ 
alytical  cantilever  beam  problem. 


The  third  group  of  problems  adds  kinetic  energy  terms 
associated  with  a  reference  frame  rotating  at  a  constant 
angular  velocity,  including  coupling  terms  (i.e.,  Corio¬ 
lis)  with  motion  relative  to  the  reference  frame.  The 
solutions  are  natural  frequencies  (as  a  function  of  rotor 
speed).  Comparisons  are  made  to  University  of  Mary¬ 
land  swept  tip  beam  experiment  data  and  an  analytical 
hinged  rotating  uniform  beam  problem. 

The  fourth  group  of  problems  adds  the  energy  terms 
associated  with  time  varying  boundary  loading.  The  so¬ 
lutions  are  periodic  response.  Comparisons  are  made 
to  a  full  scale  rotor,  flight-tested  extensively  during  the 
NASA/ Army  UH-60  Airloads  Program. 

Static  Elasticity  Problems 

Three  problems  tested  the  elastic  energy  terms  of  the  el¬ 
ement: 

1 .  The  analytical  problem  of  a  cantilever  elastica  with 
a  tip  moment. 

2.  The  experimental  “Princeton  Beam”  problem  of 
static  deflection  of  a  beam  with  strong,  nonlinear 
edgewise-flatwise-torsion  coupling. 

3.  The  analytical  problem  of  a  cantilever  elastica  with 
a  tip  load.  (The  section  describing  the  elastica  with 
tip  load  problem  depends  on  the  “Princeton  Beam” 
section  and  so  is  presented  after  it.) 

Elastica  with  Tip  Moment 

Helicopter  rotor  blades,  particularly  for  hingeless  and 
bearingless  blades  cantilevered  to  the  hub,  are  typically 
quite  flexible  and  experience  moderately  large  elastic  de¬ 
formations  for  which  conventional  linear  beam  theory  is 
often  inadequate. 

The  nonlinear  beam  element  formulation  is  based  on 
the  moderate  deformation  theory  of  Hodges  and  Dow¬ 
ell,  [13],  but  the  addition  of  a  reference  frame  attached  to 
the  root  of  each  beam  element  and  constrained  to  move 
with  the  tip  node  of  the  parent  beam  element  enables  arbi¬ 
trarily  large  deformation  to  be  accommodated,  provided 
a  sufficient  number  of  nonlinear  beam  elements  is  em¬ 
ployed.  The  elastica  problem  comparison  with  RCAS 
shows  this  approach  to  be  very  effective,  even  for  beam 
deflections  far  beyond  those  likely  to  be  encountered  in 
rotorcraft  applications.  The  elastica  problem,  described 
in  many  textbooks  (e.g.,  [14])  epitomizes  large  deflec¬ 
tions  of  beams  in  two  dimensions.  Elastica  is  the  name 
for  the  shape  of  elastic  curve  found  from  solving  the  dif¬ 
ferential  equation  (Eq.  1)  of  the  deflection  curve  when  the 


slopes  and  deflections  become  large. 
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Here  M  is  the  bending  moment,  E  is  the  modulus  of  elas¬ 
ticity,  I  is  the  moment  of  inertia,  0  is  the  angle  of  rotation 
of  the  deflection  curve,  and  s  is  the  distance  along  the 
curve. 

If  the  curvature,  K,  is  constant,  it  is 
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With  0  =  0  at  s  =  0,  the  solution  is 
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Here  p  is  the  radius  of  curvature.  (The  theory  assumes 
the  beam  is  thin  (its  thickness  h  «  p).) 

Selecting  a  coordinate  system,  the  Cartesian  coordi¬ 
nates,  rx,ry  of  any  point  on  the  deflection  curve  can  be 
calculated  as 

s 

rx=  psin0  =psin—  (4) 

P 

j1 

ry  =  p  (1  —  cos  0)  =p(l— cos-)  (5) 
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The  RCAS  model  developed  for  comparison  with  this 
closed  form  solution  was  a  uniform  cantilever  beam  with 
a  moment  applied  at  the  tip.  The  model  had  N  =  20  ele¬ 
ments,  each  of  length  /  =  1  ft,  and  an  El  =  9000 lb/ft2. 

Selecting  a  value  of  M  =  500ft-lb  resulted  in  p  = 
18  ft.  Figure  1  compares  the  closed  form  solution 
(Eqs.  4  and  5)  and  the  solution  computed  with  RCAS 
for  s  =  Oft, . . .  ,20ft.  Note  that  the  maximum  angle  0  = 
1.1 11  rad  =  .353771  =  63.66°. 

The  results  are  identical  at  the  figure’s  resolution.  Note 
that  this  is  for  a  beam  with  a  tip  deflection  of  about  56% 
L  and  a  tip  slope  of  about  64°,  values  far  beyond  the  de¬ 
flections  normally  encountered  by  a  rotor  blade. 

The  largest  discrepancies  in  Fig.  lare  of  the  order 
6  x  10-4.  Part  of  the  error  derives  from  the  numerical 
methods  used  to  solve  the  equations.  Another  portion  de¬ 
rives  from  the  element’s  use  of  polynomials  in  contrast 
to  the  trigonometric  functions  in  the  exact  solution.  This 
second  form  of  error,  which  is  intrinsic  to  the  element, 
can  be  analyzed,  giving  insight  to  the  effect  of  mesh  size, 
A,  on  accuracy. 

For  any  single  element  and  reasonable  mesh  size,  s  <C 
p.  (In  the  above  example,  s  =  L/N  =  1  ft  <C  p  =  18ft.) 
Expanding  the  sin  and  cos  terms  of  Eqs.  4  and  5  in  a 


Figure  1:  Elastica:  analytical  vs.  RCAS  for  radius  of 
curvature  p  =  18  ft  and  L  =  20  ft  with  M  =  500ft-lb, 
El  =  90001b-ft2,  and  N  =  20. 


series  gives 
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The  expressions  used  by  the  element  differ  slightly.  The 
element’s  rx  term  (which  derives  from  integration  of  a 
square  root  accounting  for  geometric  effects)  matches  the 
first  two  terms  in  Eq.  6  but  differs  from  the  third  by  being 
three  times  as  large  and  opposite  in  sign.  The  element’s  ry 
term  is  the  same  as  the  first  term  in  Eq.  7.  (This  problem 
has  no  axial  or  shear  loads  so  the  linear  term  in  rx  and 
the  cubic  term  in  ry  don’t  appear.)  The  normalized  (by 
the  analytical  solution)  difference  between  the  element 
and  analytical  solutions  at  the  tip  of  the  element  can  be 
written  (at  least  for  small  0) 

P  ~  -1  -V4  I  ,  _ 

£V  ~  30-jfr\s=L/N — 

£y  ~  T2p\s=L/N  = 


30p4(V4 
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For  the  above  example,  L/N  =  1ft  and  p  =  18  ft,  the 
error  predictions  are  about  £A  =  —2  x  10-8  and  £v  = 
3  x  10-4.  The  actual  numerical  differences  for  the  RCAS 
implementation  are  £A  =  1  x  10-6  and  £v  =  —2  x  10-3 


0  1000  2000  3000  4000  5000  6000  7000  8000  9000 

Moment  M  [ft-lb] 


Figure  2:  Elastica:  analytical  vs.  RCAS  for  radius  of  Figure  3:  Relative  intrinsic  element  error:  predicted  vs. 
curvature  p  =  3.6ft  and  L  =  20ft  with  M  =  2500 ft-lb,  actual  RCAS  error  at  the  end  of  the  first  element  with 
El  =  90001b-ft2,  and  N  =  20.  L  =  20ft,  ,v  =  1  ft,  El  =  90001b-ft2,  and  N  =  20. 


—  quite  different.  Presumably,  the  numerical  errors  as¬ 
sociated  with  the  solution  are  larger  than  the  error  asso¬ 
ciated  with  the  element  shape  functions.  This  is  for  a  sin¬ 
gle  element,  the  nonlinear  accumulation  of  errors  across 
N  elements  is  more  complex. 

The  error  estimates  (Eqs.  8  and  9)  suggest  the  effect 
of  increasing  radius  of  curvature  is  the  same  as  the  ef¬ 
fect  of  increasing  number  of  elements.  So  the  effect  of 
one  fifth  as  many  elements  on  intrinsic  error  is  the  same 
as  that  of  quintupling  the  moment.  Selecting  a  moment 
of  M  =  2500 ft-lb  resulted  in  p  =  3.6ft.  Figure  2  com¬ 
pares  the  closed  form  solution  and  the  solution  computed 
with  RCAS.  Note  that  in  this  case  the  maximum  angle 
0  =  5. 556 rad  =  1 .76871  =  318.31°,  which  is  nearly  a  full 
revolution. 

For  the  above  example,  L/N  =  1ft  and  p  =  3.6ft, 
the  error  estimates  are  £x  =  —1.98  x  10-4  and  £v  = 
6.43  x  10-3.  The  actual  (RCAS)  numerical  differences 
are  £x  =  — 1.94  x  10-4  and  ev  =  6.07  x  10-3  —  as  pre¬ 
dicted.  ry  decreases  the  radius  of  the  elastica.  Combined 
with  the  more  accurate  prediction  of  length  this  results 
in  an  over  prediction  of  0  accumulating  to  0  =  322.86°, 
an  extra  4.5°  at  the  tip.  Surprisingly,  the  0.6%  error  in 
ry  for  a  single  element  when  accumulated  over  N  =  20 
elements  leads  to  a  only  1.43%  error  in  0. 

This  error  comparison,  the  predicted  intrinsic  error  due 
to  element  shape  functions  vs.  the  actual  error  in  the 


RCAS  run,  is  presented  for  a  variety  of  applied  moments 
in  Fig.  3  —  noting  that  relative  error  shown  is: 


Actual  ry  —  Predicted  ry 
Predicted  ry 


(10) 


Note  that  if  the  total  error  includes  element  intrinsic 
error  and  solution  error,  dividing  by  the  deflection,  ry, 
amplifies  the  solution  error  for  small  moments  (and  de¬ 
flections),  explaining  the  larger  discrepancies  for  small 
moments. 

The  actual  error  is  consistent  with  the  expected  intrin¬ 
sic  error  associated  with  the  shape  functions  used.  In¬ 
terestingly,  the  error  is  less  than  9%  at  M  =  9000  ft-lb 
corresponding  to  a  1  rad  rotation  in  the  first  element  and 
over  three  complete  revolutions  for  the  N  =  20  elements. 


“Princeton  Beam”  Static  Deformations 

The  Princeton  Beam  Problem  directly  addresses  the  im¬ 
portant  nonlinear  torsion  deflections  that  accompany  the 
combined  edgewise  and  flatwise  bending  of  a  beam.  This 
torsion-bending  coupling  has  important  aeroelastic  con¬ 
sequences  for  blade  loads  and  stability.  Comparisons  be¬ 
tween  static  deformation  data  from  the  “Princeton  Beam” 
experiment  and  RCAS  are  presented  in  the  next  four 
subsections,  which  describe  the  experiment,  describe  the 


RCAS  model  used,  compare  horizontal  and  vertical  de¬ 
flections,  and  compare  twist  deflections. 

Background 

For  many  years,  a  standard  of  comparison  for  nonlinear 
beam  static  and  dynamic  response  has  been  the  “Prince¬ 
ton  Beam”  experiment.  The  data  provide  nonlinear  static 
deformation  and  natural  frequencies  for  a  long  slender 
cantilever  beam  with  a  variety  of  tip  weights.  Rotating 
the  beam  to  various  angles  induced  strong  nonlinear  cou¬ 
pling  of  edgewise-flatwise-torsion  motion.  These  exper¬ 
iments  were  performed  by  E.  H.  Dowell  and  J.  J.  Tray- 
bar  circa  1974-1975  at  the  Department  of  Aerospace  and 
Mechanical  Sciences,  Princeton  University. 

Static  data  from  the  experiment  has  been  widely  used 
(e.g.,  [15, 16]),  to  evaluate  nonlinear  beam  theories  —  not 
always  with  satisfactory  results.  Note  that  the  first  com¬ 
parisons  of  the  Hodges-Dowell  theory  with  the  experi¬ 
ments  [15]  showed  significant  deficiencies  of  the  mod¬ 
erate  deformation  nonlinear  beam  theory  at  the  largest 
experimental  loading  conditions.  Rosen  and  Fried¬ 
mann  [17],  using  an  improved,  3rd-order  nonlinear  beam 
formulation  showed  considerable  improvement,  but  this 
approach  is  still  subject  to  the  limitations  of  a  theory  de¬ 
veloped  using  an  ordering  scheme.  The  experiment  was 
carried  out  in  two  stages: 

1.  Dowell  and  Traybar  describe  the  first  stage  in  [18], 
They  measured  vertical  and  lateral  static  deforma¬ 
tion  of  three  horizontal,  cantilever,  slender,  rectan¬ 
gular  beams,  to  which  they  connected  a  variety  of 
tip  masses,  while  they  oriented  the  beams  at  a  va¬ 
riety  of  angles.  (Due  to  the  rotor  blade  focus  of 
the  experiment,  “flap”  was  used  to  refer  to  flatwise 
motion,  “lag”  or  “chordwise”  was  used  to  refer  to 
edgewise  motion,  “pitch”  was  used  to  refer  to  axial 
rotation,  and  “span”  or  “radius”  was  used  to  refer  to 
axial  length.)  They  also  measured  the  first  flapwise 
bending  and  first  lag  wise  bending  frequencies. 

2.  Dowell  and  Traybar  describe  the  second  stage 
in  [19].  From  the  beams  used  in  the  first  experiment, 
they  selected  the  beam  designated  #2  for  more  accu¬ 
rate  measurements  of  static  deformation.  They  im¬ 
proved  the  accuracy  of  their  vertical  measurements 
by  two  orders  of  magnitude  and  their  lateral  mea¬ 
surements  by  one  order  of  magnitude.  In  addition  to 
measurements  at  the  tip,  they  added  measurements 
at  intermediate  stations.  They  made  additional  mea¬ 
surements  from  which  torsional  deformation  could 
be  estimated.  As  before,  they  made  these  measure¬ 
ments  for  a  variety  of  tip  masses  and  pitch  angles. 
They  did  not  repeat  their  frequency  measurements. 


“Princeton  Beam”  RCAS  Model 

The  more  accurate  data  from  the  second  experiment  were 
used  for  comparison.  The  beam  (designated  #  2)  was  of 
length  R  =  20  in,  of  width  b  =  0.5  in,  and  of  height  h  = 
0.125  in.  The  beam  was  machined  from  7075-T651  Alu¬ 
minum  alloy,  for  which  the  material  properties  are  [20] 
—  density,  p  =  0.101  lb/in3;  Poisson’s  ratio,  v  =  .33; 
Shear  Modulus,  G  =  3.9  x  106 lb/in2;  Young’s  Mod¬ 
ulus  (averaging  the  tension  and  compression  moduli), 
E  =  10.4  x  106lb/in2. 

The  usual  formulae  were  used  for  cross-section  prop¬ 
erties  ( A  =  bh,  I  flap  =  blv’  j  12,  Iiag  =  b^hj  12).  The  tor¬ 
sional  rigidity,  GJ ,  is  a  little  more  challenging.  The 
sometimes  erroneously  used  polar  moment  of  iner¬ 
tia,  J  polar  =  5.67/z4  applies  to  a  circular  cross-section. 
The  frequently  used  membrane  analogy  (e.g.,  [21]), 
Jmembrane  =  w/z3/3  =  1.33/z4  only  becomes  accurate  for 
larger  b/h  ratios.  For  this  width  to  height  ratio,  b/h  =  4, 
Timoshenko  and  Goodier’s  result  [22]  is  preferred,  J  = 
1.124/z4  =  2.744  x  10_4in4. 

An  RCAS  model  of  the  beam  was  developed  with 
N  =  8  elements  —  which  the  following  error  analysis 
suggests  as  appropriate.  The  maximum  reported  flap  de¬ 
formation  of  the  beam  was  W  =  6 . 1 1 8  in  with  a  tip  load  of 
P  =  3.01b  and  a  pitch  angle  of  ©  =  —45°.  (The  “Prince¬ 
ton  Beam”  pitch  angle,  0,  is  measured  relative  to  a  ver¬ 
tical  chord  —  so  ©  =  0°  corresponds  to  lagwise  loading 
and  ©  =  ±90°  to  flapwise.)  The  maximum  (root)  flap- 
wise  moment  is  M  =  3.01b  x  sin (°— 45)  x  20in  which 
corresponds  (Eq.  2)  to  a  radius  of  curvature  p  =  19.95  in. 
Referring  to  Eq.  9  the  anticipated  error  at  the  end  of  the 
first  element  is  1.31  x  10-3.  (This  estimate  is  made  rec¬ 
ognizing  that  shear  and  axial  loads,  non-constant  mo¬ 
ments  and  coupling  are  present,  but  not  accounted  for  in 
Eqs.  8  and  9.)  Accumulated  over  the  N  =  8  elements,  the 
intrinsic  error  would  be  expected  to  be  of  the  order  1%. 

The  early  comparisons  were  quite  good  —  but  a  sys¬ 
tematic  pattern  was  evident  in  the  differences  between  ex¬ 
perimental  and  RCAS  results.  The  discrepancies  tended 
to  be  linear  in  both  load  and  axial  station.  This  pat¬ 
tern  is  consistent  with  compliance  in  the  support  fixture 
(a  milling  machine  type,  precision,  indexing-chuck)  or 
in  the  root  end  of  the  beam.  Regression  analysis  indi¬ 
cated  that  introducing  a  root  flap  hinge  with  a  spring  of 
26,500ft-lb/rad  in  the  RCAS  model  would  eliminate  the 
linear  trend  in  the  discrepancies. 

“Princeton  Beam”  Deflections 

The  experimenters  measured  vertical  deflection  (report¬ 
edly  to  ±.001  in)  and  horizontal  deflection  (reportedly  to 


±.01  in).  In  addition  to  the  tip  ( r  =  R),  measurements 
were  made  at  r  =  R/4,  r  =  R/2,  and  r=3R/4.  They  then 
transformed  the  measurements  to  flap  deflection,  W ,  and 
lag  deflection,  V.  It  is  tempting  to  make  comparisons  in 
these  more  natural  coordinates,  but  that  tends  to  obfus¬ 
cate  the  errors. 

The  measurements  were  always  made  relative  to  the 
unloaded  (P  =  0)  condition.  Not  relative  to  zero  defor¬ 
mation.  In  particular,  that  condition  allows  deformation 
of  the  beam  due  to  its  own  weight.  The  running  mass 
amounts  to  m  fa  6  x  10-3  lb/in  or  about  1/8  lb.  This  can 
be  significant  for  the  smaller  loads,  especially  P  =  1  /21b. 
Accurate  comparison  requires  subtracting  the  RCAS  so¬ 
lution  for  no  tip  load  (P  =  0)  but  including  the  deforma¬ 
tion  due  to  the  beam’s  own  weight  from  the  solution  with 
the  specified  tip  load  (and  the  beam’s  own  weight).  De¬ 
flections  of  the  order  1/8  in  were  predicted  for  the  P  =  0 
condition,  corresponding  to  Z/L  =  0.006,  about  the  size 
of  the  larger  symbols  on  the  following  figures  (Figs.  4 
and  5). 

A  comparison  of  the  RCAS  and  “Princeton  Beam’’  hor¬ 
izontal  deflections  for  a  P  =  3.01b  tip  load  is  presented 
in  Fig.  4.  (Recall  the  horizontal  deflection  is  flapwise  at 
0  =  0°  and  lagwise  at  0  =  90° .)  A  comparison  of  the 
RCAS  and  “Princeton  Beam”  vertical  deflections  for  a 
P  =  3.01b  tip  load  is  presented  in  Fig.  5.  (Also  recall  the 
vertical  deflection  is  lagwise  at  0  =  0°  and  flapwise  at 
0  =  90° .) 


“Princeton  Beam”  Twist 

Although  there  are  no  applied  torsion  moments,  per 
se.  internal  torsional  moments  are  an  inherent  nonlinear 
consequence  of  the  simultaneous  edgewise  and  flatwise 
bending.  To  facilitate  twist  measurements,  the  experi¬ 
menters  fastened  Lrr  =  6  in  reference  rods  perpendicular 
to  the  beam  (in  the  edgewise  direction)  at  r  =  R/4,  R/2, 
3R/4,  and  R.  In  the  absence  of  tip  load,  the  horizontal 
and  vertical  distances  between  the  tips  of  the  rods  are  la¬ 
beled  XpISTANCE  and  ZpISTANCE.  In  the  presence  of  tip 
load,  and  the  resulting  elastic  twist,  the  same  measure¬ 
ments  are  labeled  X'DISTANCE,  and  Z'DISTANCE.  With  or 
without  load,  they  are  nominally: 

X, distance  «  6  sin  ©inches  (11) 

Zdistance  «  6 cos  ©inches  (12) 

The  experimenters  presented  twist  calculated  in  two 
ways  —  the  arc  tangent  method  (better  when  0  «  0°) 
and  the  arc  cosine  method  (better  when  0  m  90°).  The 


0  10  20  30  40  50  60  70  80  90 

"Pitch"  Angle  0  [deg] 


Figure  4:  “Princeton  Beam”  horizontal  deflections,  X: 
experiment  vs.  RCAS  for  P  =  3.01b  and  N  =  8  with 
R  =  20in,  b  =  0.5 in,  and  h  =  0.125 in. 
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Figure  5:  “Princeton  Beam”  vertical  deflections,  Z:  ex¬ 
periment  vs.  RCAS  for  P  =  3.01b  and  N  =  8  with  R  = 
20in,  b  =  0.5 in,  and  h  =  0.125 in. 


Figure  6:  “Princeton  Beam”  twist  angle  (arc  tangent 
method),  <j>:  experiment  vs.  RCAS  for  P  =  3.01b  and 
N  =  8  with  R  =  20in,  b  =  0.5 in,  and  li  =  0.125 in. 

formulae  they  presented  for  the  two  methods  are: 

$  =  <>  -  <>0  =  arctan  X™tance  _  arctan  ^distance 

2 DISTANCE  2 DISTANCE 

(13) 

$  =  9  -  4>0  =  arccos  Distance  _  arccos  ^distance 
6.0  6.0 

(14) 

The  arc  tangent  method  suffers  from  a  singularity  at 

2 DISTANCE  =0or  ^DISTANCE  =  °>  both  °f  which  WOllld 

be  expected  to  occur  at  ©  =  90°. 

Rigid  rods  were  included  in  the  RCAS  model  equiv¬ 
alent  to  those  in  the  experiment.  The  calculations  in 
Eqs.  13  and  14  were  carried  out  using  an  RCAS  calcu¬ 
lation  for  an  unloaded  beam,  P  =  0,  (but,  as  noted  above 
(“Princeton  Beam”  Deflections  section)  not  zero  defor¬ 
mation  due  to  running  mass)  to  calculate  Distance  and 
2 distance ■  A  comparison  of  the  results  obtained  with 
the  arc  tangent  method,  Eq.  13,  with  the  reported  results 
is  presented  in  Fig.  6.  Similarly,  a  comparison  of  the  re¬ 
sults  obtained  with  the  arc  cosine  method,  Eq.  14,  with 
the  reported  results  is  presented  in  Fig.  7. 

One  is  immediately  struck  by  the  apparently  spurious 
results  at  ©  =  0  using  the  arc  cosine  method  (Fig.  7). 
The  values  are  actually  a  correct  result  from  application 
of  Eq.  14.  The  projection  of  the  reference  rod  in  the  Z 
direction  can  be  reduced  by  either  a  pitch  rotation,  ©, 
or  a  rotation  about  the  X  axis  «  dZ/dr.  At  ©  =  0,  the 


Figure  7:  “Princeton  Beam”  twist  angle  (arc  cosine 
method),  <j>:  experiment  vs.  RCAS  for  P  =  3.01b  and 
N  =  8  with  R  =  20in,  b  =  0.5  in,  and  h  =  0.125  in. 

foreshortening  is  a  result  of  beam  axis  bending  —  not 
pitch  and  twist.  The  slope  of  the  beam  axis  is  plotted  on 
the  same  figure  (Fig.  7). 

Given  the  scatter  in  the  experimental  data,  it  may  be 
instructive  to  estimate  the  effect  of  measurement  error 
on  the  calculated  twist,  <j).  Assuming  the  beam  is  ori¬ 
ented  at  a  twist  angle,  <f>  w  ©  or  <j)(J  «  ©,  and  that  there  are 
unknown  measurement  errors,  Ey  and  £j  or  evo  and  e_o, 
equations  Eqs.  13  and  14,  with  Eqs.  11  and  12  give: 


r  6sin(b-(-£y  6sin(bo+eTo 

ffl  «  arctan - arctan - — 

ocoscp  +  e^/  6cos(|)o  +  £zo 

(15) 

~  Gcostb  +  e^  6cos(|)o+£7o 

o  «  arccos - 1 —  arccos - - - — 

v  6.0  6.0 

(16) 

Noting  the  derivatives: 

d  6sind)+x  cosd) 

arctan  =  , 

ox  6  cos  (p  6 

(17) 

d  6  costb  -\-x  —1 

—  arccos - = - - 

ox  6.0  6sin(p 

(18) 

Neglecting  £,  in  comparison  to  £x,  the  errors  in  the  two 
calculations  are: 

cos©  cos© 

~  g  ey  6  e-v° 

(19) 

^~6sin©£;:'  6sin©^° 

(20) 

For  ©  =  15°,  30°,  and  45°,  the  sensitivities  are  about 
9°  /in,  8°  /in,  7° /in,  respectively,  for  the  arc  tangent 
method;  and  37° /in,  19° /in,  14° /in,  respectively,  for  the 
arc  cosine  method.  The  experimenters  reporting  mea¬ 
surement  accuracies  of  ±.01  in  laterally,  and  ±.001  in 
vertically.  The  statistical  combination  of  the  errors  in 
the  eight  measurements  (two  for  XEISTA NCE,  ■  ■  ■ )  multi¬ 
plied  by  the  sensitivities  yield  about  ±0.18°,  ±0.16°,  and 
±0.14°  for  the  arc  tangent  method  and  ±0.07°,  ±0.04°, 
and  ±0.03°  for  the  arc  cosine  method  at  ©  =  15°,  30°, 
and  45°,  respectively. 

To  limit  the  clutter,  these  error  estimates  are  illustrated 
in  Figs.  6  and  7  for  only  station  r/R  =  0.25.  For  each 
angle,  ©,  the  size  of  the  error  bar  is  the  same  at  all  sta¬ 
tions.  The  data  scatter  using  the  arc  tangent  method  is 
consistent  with  the  predicted  error  estimates  and  the  re¬ 
ported  measurement  error;  the  scatter  using  the  arc  cosine 
method  is  less  so. 

Since  the  reported  calculations  seem  to  obscure  the 
comparison,  the  calculation  was  inverted.  Assuming 
Eqs.  11  and  12  are  exact  (instead  of  approximate  when 
X distance  is  replaced  with  XEISTANCE  (and  Zd1STance 
with  Z/)ISTANCE),  Eqs.  13  and  14  can  be  solved  for 
X distance  an£l  Z 'distance  ■  Calculating  the  difference  be¬ 
tween  the  measurement  and  the  nominal  values  (Eqs.  1 1 
and  12)  and  normalizing  by  the  length  of  the  reference 
rod,  Lrr  =  6  in,  the  twist  related  measurement  is  com¬ 
pared  with  the  same  calculation  using  RCAS  results  in 
Figs.  8  and  9. 

Again,  to  limit  the  clutter,  these  error  estimates  are  il¬ 
lustrated  in  Figs.  8  and  9  for  only  station  r/R  =  0.25.  For 
each  angle,  ©,  the  size  of  the  error  bar  is  the  same  at  all 
stations.  The  “Princeton  Beam”  data,  the  error  bounds 
and  the  RCAS  results  are  all  reasonably  consistent. 

Elastica  with  Tip  Load 

Another  problem  for  which  an  analytical  solution  exists 
is  the  tip  loaded  elastica.  The  solution  is  in  the  form  of  a 
transcendental  equation  involving  elliptic  functions.  So¬ 
lutions  are  available  from  a  variety  of  sources  —  those 
from  reference  [14]  are  used  for  comparison  here.  Nor¬ 
malized  elastica  angle  of  rotation,  0/,/lj,  vertical  deflec¬ 
tion,  8 V/L,  and  horizontal  deflection,  8/,  / L  (all  at  the  tip) 
vs.  normalized  loads,  PL1  /El,  are  tabulated  op.  cit.  for 
a  cantilever  beam  of  length  L,  modulus  of  elasticity  E, 
and  moment  of  inertia  /,  subjected  to  a  (downward)  tip 
load  P.  Note  that  horizontal  deflection,  8/,  for  this  prob¬ 
lem  refers  to  the  shortening  of  the  projection  of  the  beam 
axis  in  the  axial  direction,  not  lateral  deflection  as  with 
the  “Princeton  Beam”.) 

The  RCAS  model  developed  for  the  “Princeton  Beam” 
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Figure  8:  “Princeton  Beam”  twist:  normalized  relative 
lateral  deflection,  X' /Lrr  —  sin®:  experiment  vs.  RCAS 
for  P  =  3.01b  and  N  =  8  with  R  =  20 in,  b  =  0.5 in,  and 
h  =  0.125  in. 
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Figure  9:  “Princeton  Beam”  twist:  normalized  relative 
vertical  deflection,  cos®  —  Z' /Lrr:  experiment  vs.  RCAS 
for  P  =  3.01b  and  N  =  8  with  R  =  20 in,  b  =  0.5 in,  and 
/;  =  0.125  in. 


Figure  10:  Tip  Loaded  Elastica:  Rotation  Angle,  0 j,  vs. 
RCAS  for  L  =  20in,  N  =  8,  P  =  0.5...  101b  and  El  = 
8.463  x  102  or  1.354  x  104 lb/in2. 

comparison,  discussed  above,  was  reused  since  it  allowed 
direct  comparison  to  the  experimental  results.  The  flap 
hinge  and  spring  (which  modeled  root  compliance)  was 
the  only  feature  removed.  In  particular,  shear  and  axial 
deflections  were  allowed  but  were  expected  to  be  negli¬ 
gible.  Running  mass  was  retained,  although,  as  discussed 
previously  (“Princeton  Beam”  Deflections  section)  the 
unloaded  solution  (P  =  0)  was  subtracted  from  the  RCAS 
solutions.  This  model  only  approximated  the  analytical 
problem,  so  exact  comparison  was  not  anticipated. 

Restricting  this  N  =  8  element,  L  =  20  in  model  to 
0  =  0°  and  90°  eliminated  the  edgewise-flatwise-torsion 
coupling.  At  0°,  the  bending  is  lagwise  with  El  = 
1.354  x  104 lb/in2  and  at  90°,  the  bending  is  flapwise 
with  El  =  8.464  x  102  lb/in2. 

The  RCAS  model  was  run  for  tip  loads  of  P  = 
0.5, . . . ,  101b  corresponding  to  normalized  loads  of  0.24- 
4.7  flapwise  and  0.015-0.30  edgewise.  A  comparison 
of  analytical  and  RCAS  results  is  presented  in  Fig.  10 
for  normalized  tip  bending  rotation.  Fig.  1 1  for  normal¬ 
ized  vertical  displacement,  and  Fig.  12  for  normalized 
axial  displacement.  At  the  largest  load,  P  =  101b,  the 
beam  tip  angle  is  0*  «  68°,  the  vertical  displacement 
is  8V  m  14in  =  ,7L,  and  the  horizontal  displacement  is 
8V  «  7.6 in  =  .38L. 

The  RCAS  edgewise  data  points  only  appear  as  a 
thickening  of  the  line  in  the  lower  left  corner  of  each 


Figure  1 1 :  Tip  Loaded  Elastica:  Vertical  Displacement, 
8V  vs.  RCAS  for  L  =  20in,  N  =  8.  P  =  0.5...  101b  and 
El  =  8.463  x  102  or  1.354  x  104 lb/in2. 


Figure  12:  Tip  Loaded  Elastica:  Horizontal  Displace¬ 
ment,  8/,  vs.  RCAS  for  L  =  20in,  N  =  8.  P  =  0.5 ...  101b 
and  El  =  8.463  x  102  or  1.354  x  104 lb/in2. 


Figure  13:  Tip  Loaded  Elastica:  Vertical  Displacement, 
8V  vs.  RCAS  and  vs.  “Princeton  Beam”  for  L  =  20 in, 
N  =  8,  P  =  0.5  ...  101b  and  El  =  8.463  x  102  or  1 .354  x 
104  lb/in2. 

graph.  The  data  from  the  “Princeton  Beam”  experiment 
are  superimposed  on  Fig.  1 1  and  presented  at  a  magnified 
scale  in  Fig.  13. 

The  analytical,  “Princeton  Beam”,  and  RCAS  results 
are  quite  consistent,  confirming  the  expectations  that 
shear  and  axial  deformation  were  negligible.  It  also  sup¬ 
ports  the  expectation  that 


presence  of  the  substantial  edgewise-flatwise-torsion 
static  deformations  discussed  above. 

For  the  static  measurements,  the  experimenters  sus¬ 
pended  a  weight  pan  by  a  thread  connected  to  a  screw 
in  the  end  of  the  beam  and  added  weights  to  produce  the 
desired  static  load.  They  had  to  use  a  different  approach 
for  dynamics.  Judging  from  a  photo  and  a  drawing,  they 
machined  a  hollow  cylindrical  mass  with  a  central  web, 
which  was  slid  over  the  end  of  the  beam  and  fastened, 
through  the  web,  to  the  end  of  the  beam.  The  mass  shown 
was  of  length  about  3/4  of  the  diameter.  Other  than  the 
photo  and  drawing,  no  information  was  available  regard¬ 
ing  the  dimensions  of  the  weights. 

To  provide  a  rough  estimate  of  rotary  inertia,  the  in¬ 
ertia  was  assumed  to  be  that  of  a  hollow  cylinder  (Ixx  = 
4,  =  m(3R2  +  3 r2  +  /z2) / 12)  of  the  specified  weight  and 
of  the  same  7075-T651  Aluminum  alloy.  The  mass  and 
volume  are  thus  m  =  P/g  =  pV.  The  length  was  assumed 
to  be  h  =  1.5\/3V/(1.57t).  The  inner  radius  was  taken 
to  be  r  =  ,025  ft.  The  outer  radius  was  then  calculated 
as  \fv /  (nil)  +  r2.  The  result  is  a  hollow  cylinder,  of  the 
correct  weight,  large  enough  to  fit  over  the  beam  with  a 
little  clearance,  and  whose  length  (h  «  1 .5 R)  is  within  a 
few  percent  of  3/4  of  the  outer  diameter. 

A  comparison  of  “Princeton  beam”  experimental  and 
RCAS  results  is  presented  in  Fig.  14  for  tip  masses 
weighing  P  =  21b  and  P  =  31b.  Figure  15  presents  the 
flapwise  results  at  a  different  scale.  The  results  show 
the  effect  of  static  equilibrium  deflections  of  the  beam 
on  the  natural  frequencies  as  well  as  the  nonlinear  cou¬ 
pling  due  to  the  combined  bending  and  torsion  for  the 
non-zero  pitch  angles.  The  quite  acceptable  agreement 
validates  the  theoretical  basis  of  RCAS. 


Fixed  Frame  Dynamics  Problems  Cantilever  Beam  Frequency  Comparisons 


Two  problems  tested  the  kinetic  energy  terms  of  the  ele¬ 
ment  for  motion  relative  to  a  fixed  reference  frame: 

1 .  The  experimental  “Princeton  Beam”  problem  of  fre¬ 
quencies  at  static  equilibrium  of  a  beam  with  strong, 
nonlinear  edgewise-flatwise-torsion  coupling. 

2.  The  analytical  problem  of  the  eigenvalues  of  a  can¬ 
tilever  beam. 


In  the  absence  of  a  tip  mass,  the  results  from  the  pre¬ 
vious  (“Princeton  Beam”  Frequencies  section)  admit  an 
analytical  solution.  Many  textbooks  (e.g.,  [23])  provide 
the  natural  frequencies  for  a  uniform  cantilever  beam  as: 


C0„ 


(21) 


“Princeton  Beam”  Frequencies 

As  previously  noted  (“Princeton  Beam”  Static 
Deformations  section),  frequency  data  was  gathered 
during  the  “Princeton  Beam”  experiment,  although 
previous  investigators  have  made  very  little  use  of  this 
data.  These  measurements  provide  the  first  flapwise 
and  chordwise  natural  frequencies  for  vibrations  in  the 


where  ([3/)i  ~  1.875  is  the  reported  solution  to  the  tran¬ 
scendental  equation  for  the  first  natural  frequency. 

The  analytical  solution  (using  the  RCAS  value  for 
gravity,  g  =  32.17405  ft/sec2  to  calculate  mass),  the 
“Princeton  Beam”  experimental  frequencies,  and  the  fre¬ 
quencies  RCAS  calculated  for  the  same  conditions  (0  = 
0°  or  90°  and  P  =  0  are  tabulated  for  the  first  flatwise  and 
first  edgewise  frequencies  in  Table  1 . 
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Figure  14:  “Princeton  Beam”  natural  frequencies  of  first 
flapwise  and  chordwise  modes:  experiment  vs.  RCAS  for 
P  =  2 lb  and  P  =  3 lb  and  N  =  8  with  R  =  20in,  b  =  0.5  in, 
and  h  =  0.125  in. 
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Figure  15:  “Princeton  Beam”  natural  frequencies  of  first 
flapwise  mode:  experiment  vs.  RCAS  for  P  =  21b  and 
P  =  3  Ih  and  N  =  8  with  R  =  20  in,  b  =  0.5 in,  and  li  = 
0.125  in. 


Cantilever  Beam  Natural  Frequencies 


Model 

Flatwise 

Edgewise 

Analytical 

10.0642 

40.2570 

RCAS  0  =  0° 

10.0605 

40.0719 

RCAS  0  =  90° 

10.0602 

40.0748 

10.154 

Princeton  0  =  0° 

10.152 

10.143 

41.143 

Table  1:  Cantilever  beam  natural  frequencies:  analytical 
vs.  RCAS  vs.  “Princeton  Beam”  for  A  =  8  with  R  =  20in, 
b  =  0.5  in,  and  h  =  0.125  in. 

Rotating  Frame  Dynamics  Problems 

Two  problems  tested  the  kinetic  energy  terms  of  the  ele¬ 
ment  for  motion  relative  to  a  reference  frame  rotating  at 
a  constant  angular  velocity: 

1 .  The  experimental  frequencies  of  a  swept  tip,  rotat¬ 
ing  beam,  measured  in  the  University  of  Maryland 
vacuum  chamber. 

2.  The  analytical  problem  of  the  natural  frequencies  of 
a  hinged  rotating  beam. 

University  of  Maryland  Rotating  Beam 
Frequencies 

Comparisons  between  natural  frequency  data  from  a  Uni¬ 
versity  of  Maryland  swept  tip  beam  experiment  and 
RCAS  are  presented  in  the  next  three  subsections.  The 
experiment  is  described,  the  RCAS  model  is  described 
and  comparisons  of  natural  frequencies  are  presented. 

Swept  Tip  Background 

The  effect  of  blade  tip  sweep  angle  on  rotor  blade  struc¬ 
tural  dynamics  was  a  nontrivial  problem  for  rotorcraft 
when  such  blades  were  introduced  for  aerodynamic  and 
aeroelastic  reasons.  Celi  and  Friedmann  [24]  provided 
one  of  the  first  satisfactory  treatments  of  the  problem  al¬ 
though  good  experimental  data  for  fundamental  valida¬ 
tion  swept  tip  blade  frequencies  was  scarce  until  rotating 
blade  data  was  obtained  in  the  University  of  Maryland 
vacuum  chamber.  This  set  of  experimental  data  was  ob¬ 
tained  by  the  Department  of  Aerospace  Engineering  at 
the  University  of  Maryland.  The  “University  of  Mary¬ 
land”  data  [25]  provides  in  vacuo ,  rotating  frequencies  of 
a  blade  with  the  outboard  16%  swept  at  various  angles. 
Although  the  experiment  included  composite  as  well  as 
aluminum  blades,  only  the  results  obtained  with  the  alu¬ 
minum  blade  are  used  for  comparison. 


The  geometry  is  slightly  more  complex  than  the  prob¬ 
lems  discussed  above.  The  experimenters  used  uniform, 
untwisted  rotating  beams  with  a  16%,  0°  to  45°  swept 
tip.  They  were  spin  tested  in  the  University  of  Maryland 
vacuum  chamber. 

Experimental  frequency  measurements  have  been 
compared  with  a  finite  element  analysis  by  Epps  and 
Chandra  [25]  as  well  as  with  a  more  accurate  theory  by 
Hodges,  Shang,  and  Cesnick  [26],  Although  these  the¬ 
ories  were  reasonably  successful,  neither  study  included 
static  gravity  deflections  and  there  has  been  some  confu¬ 
sion  about  the  modal  content  of  the  analytical  and  exper¬ 
imental  results.  In  both  papers,  [25,  26],  frequencies  are 
plotted  as  a  function  of  sweep  angle  at  the  test  angular 
velocities. 

The  experimental  data  used  for  the  comparison  de¬ 
scribed  below  was  scaled  from  the  figures  in  [25].  The 
data  was  for  what  was  labeled  the  first  five  flap  modes 
and  the  first  torsion  mode.  The  rotor  speeds  presented 
were  for  0,  500,  and  750  RPM.  The  tested  beams  had 
sweep  angles  of  0°,  15°,  30°,  and  45°. 

Swept  Tip  RCAS  Model 

The  unswept  rotor  radius  was  R  =  40  in,  which  included 
a  hub  of  radius  r/,,^  =  2.5  in  =  .0625 /?,  a  beam  of  length 
rswp  —  rim},  =  31. 5  in  =  .7875 /?,  and  a  tip  section  of  length 
rtjp  —  rswp  =  6  in  =  .15/?.  Both  the  beam  section  and 
the  tip  section  were  of  width  b  =  1.0  in,  and  of  height 
h  =  0.063  in.  The  aluminum  alloy  is  not  specified  in  the 
paper,  so  it  was  assumed  to  be  7075  with  the  previously 
noted  (“Princeton  Beam”  RCAS  Model  section)  material 
properties  and  cross-section  formulae. 

Again,  torsional  rigidity  is  more  challenging.  This 
width  to  height  ratio,  b/a  =  1/0.063  =  15.873,  is  not 
tabulated  in  references  such  as  Timoshenko  and  Good- 
ier  [22] .  However,  if  the  tabulated  values  are  plotted  with 
an  reciprocal  scale  (Fig.  16),  a  value  of  Ki  =  .320  can 
be  inferred.  The  corresponding  /  =  K\lr'h  =  8.0015  x 
10-5in4. 

An  RCAS  model  of  the  beam  was  developed  with 
N  =  14  elements  to  represent  the  beam  section  and  N  =  3 
elements  to  represent  the  tip. 

Swept  Tip  Frequency  Comparisons 

Each  figure  in  the  papers  [25,  26]  presented  experimental 
data  and  theoretical  curves  for  a  single  mode’s  frequency 
as  a  function  of  sweep  angle  at  the  three  rotor  speeds. 
The  experimenters  presented  six  figures  which  they  la¬ 
beled  flap  modes  1-5  and  torsion.  All  experimental  data 
presented  here  is  labeled  according  to  their  definitions. 


Figure  16:  Reciprocal  Scale  Plot  of  Constants  for  Torsion 
of  a  Rectangular  bar  from  Timoshenko  and  Goodier  [22] . 

The  interpretation  of  the  modal  frequency  variations 
with  blade  tip  sweep  angle.  A,  and  with  rotor  speed, 
G,  are  facilitated  by  comparing  frequency  vs.  sweep  an¬ 
gle  and  frequency  vs.  rotor  speed  fan  plots.  A  com¬ 
parison  of  non-rotating  (G  =  0)  RCAS  frequencies  vs. 
experimental  data,  without  gravity  is  presented  in  Fig.  17. 
A  comparison  of  unswept  (A  =  0)  RCAS  frequencies  vs. 
experimental  data,  without  gravity  is  presented  in  Fig.  18. 

The  first  eight  flap  modes  are  easily  identified.  The  first 
torsion  and  the  first  two  lag  modes  are  also  clear,  showing 
very  little  increase  with  frequency  due  to  centrifugal  stiff¬ 
ening.  The  RCAS  frequencies  are  identified  as  f  1— f4  and 
f8  for  the  first  four  and  eighth  flap  modes;  11  and  12  for 
the  first  two  lag  modes;  cl-c3  for  three  coupled  modes 
which  at  A  =  0  and  G  =  0  are  the  fifth  through  seventh 
flap  modes;  and  c4  for  the  fourth  coupled  mode,  which  at 
A  =  0  and  G  =  0  is  the  first  torsion  mode. 

When  gravity  is  included  in  the  RCAS  calculations, 
the  non-rotating  and  unswept  frequencies  are  shown  in 
Figs.  19  and  20.  The  effects  of  gravity  are  significant 
only  at  the  lowest  rotor  speeds,  as  would  be  expected 
with  a  centripetal  acceleration  of  nearly  640  g  at  the  tip 
for  G„„m  dominating  the  acceleration  of  gravity  at  other 
rotor  speeds.  At  the  lowest  rotor  speeds,  the  weight  of 
the  blade  results  in  flapwise  or  primarily  flapwise  (when 
A  ^  0)  static  deflection.  This  static  deflection  has  only 
a  small  effect  on  the  blade  flap  frequencies.  However, 


Figure  17:  Swept  tip  natural  frequencies:  experiment  vs.  Figure  19:  Swept  tip  natural  frequencies:  experiment  vs 
RCAS  for  N  =11  and  fl  =  0  with  no  gravity.  RCAS  for  N  =  17  and  11  =  0  with  gravity. 


Normalized  Rotational  Speed,  £1  /  ^nom 


Figure  18:  Swept  tip  natural  frequencies:  experiment  vs.  Figure  20:  Swept  tip  natural  frequencies:  experiment  vs 
RCAS  for  N  =  17  and  A  =  0°  with  no  gravity.  RCAS  for  N  =  17  and  A  =  0°  with  gravity. 


Normalized  Rotational  Speed,  Q.  /  ^nom 

Figure  21:  Swept  tip  natural  frequencies:  experiment  vs. 
RCAS  for  N  =  17  and  A  =  15°  with  gravity. 

as  evident  in  the  “Princeton  Beam”  analysis,  the  edge¬ 
wise  frequencies  are  significantly  influenced  by  the  static 
deflection  induced  coupling  of  the  edgewise  and  torsion 
motion.  The  effects  of  gravity  were  not  included  in  the 
analytical  results  of  references  [25,  26],  however,  this  did 
not  affect  the  analytical  experimental  comparisons  since 
only  flapwise  and  torsion  mode  experimental  measure¬ 
ments  were  obtained.  As  a  result,  the  predominant  ef¬ 
fects  of  gravity  cannot  not  be  compared  with  the  present 
calculations. 

A  more  subtle  effect  of  gravity  is  the  coupling  of  tor¬ 
sion  with  the  fifth  flap  mode.  On  the  left  side  of  Fig.  20, 
cl  is  the  fifth  flap  mode  and  c4  is  the  first  torsion  mode, 
but  on  the  right  side  the  roles  are  reversed.  The  curves 
which  crossed  in  Fig.  1 8  only  closely  approach  each  other 
in  Fig.  20.  A  similar  interaction  occurs  between  the  sec¬ 
ond  lag  mode  and  the  sixth  flap  mode.  (However  the  cou¬ 
pled  lag  mode  will  still  be  labeled  12  to  facilitate  focusing 
on  flap  coupling.) 

With  the  introduction  of  a  A  =  15°  sweep  angle  the 
frequency  as  a  function  of  rotor  speed  is  presented  in 
Fig.  21.  The  first  torsion  mode  is  strongly  affected.  The 
mode  labeled  cl,  is  returning  to  the  characteristics  of  a 
fifth  flap  mode.  The  modes  labeled  c2  and  c4  are  coupled 
torsion  and  sixth  flap  modes  with  some  coupling  to  11, 
the  lag  mode. 

At  the  maximum  experimental  sweep  angle,  A  =  45°, 
(Fig.  22)  illustrates  the  frequencies  behavior  with  respect 


Figure  22:  Swept  tip  natural  frequencies:  experiment  vs. 
RCAS  for  N  =  17  and  A  =  45°  with  gravity. 

to  rotor  speed.  The  mode  c4  is  taking  on  more  of  the  char¬ 
acteristics  of  the  sixth  flap  mode  (especially  at  higher  ro¬ 
tor  speeds).  The  modes  c2  and  c3  are  coupled  first  torsion 
and  seventh  flap  modes. 

With  the  addition  of  a  plot  of  modal  frequency  vs. 
sweep  angle  at  the  nominal  rotor  speed,  Q.nom  (Fig.  23) 
in  conjunction  with  Fig.  22,  Fig.  19,  and  Fig.  20,  it  is 
possible  to  trace  each  eigenvalue  around  a  closed  cir¬ 
cuit  and  insure  consistent  identification.  The  evolution 
of  the  torsion  mode  as  it  couples  successively  with  the 
fifth,  sixth,  and  seventh  flap  modes  is  readily  apparent 
in  Fig.  23  The  experimental  frequency  measurements  re¬ 
flect  only  the  first  six  flap  modes,  that  is,  the  torsion  mode 
frequency  for  45  deg  sweep  angle  was  above  the  range 
of  the  measurements  (and  the  reported  measurement  was 
most  likely  the  sixth  flap  mode.  These  observations  differ 
in  some  respects  from  the  designation  of  modes  given  in 
references  [25,  26]. 

The  comparisons  indicate  that  RCAS  correctly  models 
the  effects  of  sweep  angle  on  the  rotating  blade  coupled 
flap  and  torsion  frequencies. 

Hinged  Rotating  Beam  Comparisons 

As  noted  in  Aeroelasticity  [27]  and  other  texts,  the  eigen- 
solutions  for  a  rotating  beam,  pinned  at  the  axis  of  ro¬ 
tation,  can  be  analytically  represented  in  series  form  us¬ 
ing  Duncan  polynomials.  Harris  [28],  by  specifying  a 


0  5  10  15  20  25  30  35  40  45 

Tip  Sweep  Angle,  A  [deg] 


Figure  23:  Swept  tip  natural  frequencies:  experiment  vs. 
RCAS  for  A  =17  and  £2  =  Q.nom  with  gravity. 


solution  with  a  finite  number  of  rational  coefficients,  de¬ 
fined  a  problem  with  that  exact  solution.  He  proposed 
that  problem  be  used  as  an  analytical  test  case  for  val¬ 
idating  finite  element  computer  programs.  The  rele¬ 
vant  parameters  for  the  eigenvalue  problem  are  R  =  28  ft, 
Q.  =  Vtip/R  =  650 /28  rad/sec,  El  =  5R3  =  109,7601b-ft2, 
and  m  =  168/845  slug/ft.  The  corresponding  dimen¬ 
sionless  reference  frequency  parameter  ([27])  is  1  /K  = 
mQrR4/EI  =  1/600.  Solution  of  the  transcendental 
equation  associated  with  the  analytical  problem,  retain¬ 
ing  the  first  70  terms  in  the  series,  has  been  carried  out 
by  the  author  of  the  problem  ([28])  and  made  available 
through  a  personal  communication  [29]. 

RCAS  models  were  developed  with  N  =  10  and  N  = 
20  elements.  Motion  was  restricted  to  flap  displacement 
and  rotation.  The  root  was  hinged  at  the  rotational  axis. 
A  comparison  of  the  RCAS  and  analytical  results  [29]  is 
presented  in  Table  2. 

As  anticipated,  errors  (assuming  the  analytical  solution 
is  exact)  using  finite  elements  are  consistently  small  and 
decrease  when  more  elements  are  used.  For  the  second 
mode,  the  error  was  0.005%  for  N  =  10  and  .001%  for 
N  =  20.  For  the  tenth  mode,  the  error  was  2%  for  N  =  10 
and  0.2%  for  N  =  20. 


Natural  Frequencies 

Rotating,  Pinned-Free  Beam 

Mode 

RCAS  N=10 

RCAS  N=20 

Analytical 

1 

0.999999 

0.999999 

1.0 

2 

2.55726 

2.55714 

2.55711 

3 

4.58062 

4.58006 

4.57997 

4 

7.24717 

7.24458 

7.24431 

5 

10.5852 

10.5742 

10.5733 

6 

14.6469 

14.6097 

14.6067 

7 

19.4912 

19.3883 

19.3798 

8 

25.1758 

24.9371 

24.9164 

9 

31.7440 

31.2763 

31.2299 

10 

39.1465 

38.4232 

38.3444 

Table  2:  Hinged  rotating  beam  natural  frequencies:  an¬ 
alytical  [29]  vs.  RCAS  for  N  =  10  and  N  =  20  element 
models. 


Periodic  Boundary  Load  Problems 

One  problems  tested  the  work  terms  associated  with  time 
varying  boundary  loading.  The  flight-test  loads  of  a  real 
world,  full  scale,  UH-60  aircraft  rotor  blade  experiencing 
periodic  loading,  during  steady  level  flight. 

UH-60  Blade  Loads 

The  prediction  of  rotor  blade  loads  in  response  to  un¬ 
steady  aerodynamic  forces  is  one  of  the  most  difficult 
problems  in  rotorcraft  analysis.  For  most  flight  con¬ 
ditions,  complex  fluid  flow  phenomena  prevent  aerody¬ 
namic  lift,  drag,  and  moments  from  being  calculated  ac¬ 
curately  enough  to  determine  the  dynamic  response  and 
structural  loads  of  the  rotor  blades  with  accuracy  needed 
for  design  purposes.  Even  if  accurate  aerodynamic  anal¬ 
yses  were  available,  the  complex  structural  dynamics  of 
a  nonuniform,  flexible,  rotating  beam  remains  a  very 
daunting  problem  with  little  evidence  to  demonstrate  that 
accurate  solutions  can  be  achieved  with  current  analy¬ 
sis  methods.  This  problem  involves  dynamic  coupling 
of  the  fuselage,  rotor,  controls,  and  engine  drive  train 
components.  The  blades  involve  large  rigid  body  pitch, 
flap,  and  lead-lag  rotations,  as  well  as  significant  elas¬ 
tic  deformations  in  flapwise,  chordwise  bending,  and  tor¬ 
sion.  The  large  rigid  body  blade  rotations  involve  nonlin¬ 
ear  kinematic  effects  and  the  blade  elastic  deflections  in¬ 
volve  nonlinear  coupling  of  bending  and  torsion  motions. 
Therefore,  the  full  problem  involves  many  complicated 
dynamic  and  structural  dynamics  effects. 

The  present  problem  affords  gaining  insight  into  the 
accuracy  of  rotor  blade  structural  dynamics  analysis  by 
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applying  experimentally  measured  airloads  to  a  rotor 
blade  model  and  comparing  the  resulting  blade  loads  with 
corresponding  experimental  blade  loads.  This  opportu¬ 
nity  is  available  as  a  result  of  the  full-scale  flight  test 
investigation  of  the  NASA/ Army  UH-60  Airloads  Pro¬ 
gram  reported  in  [30],  An  extensively  instrumented  UH- 
60  aircraft  was  used  to  comprehensively  measure,  among 
other  things,  blade  motion,  surface  pressures,  and  struc¬ 
tural  blade  loads  for  a  variety  of  flight  conditions.  The 
present  problem  describes  the  initial  results  of  an  investi¬ 
gation  to  apply  the  RCAS  computer  program  to  calculate 
rotor  blade  loads  in  response  to  the  measured  airloads  and 
in  turn  compare  the  calculated  blade  loads  with  the  corre¬ 
sponding  measured  loads.  In  this  way  the  accuracy  of  the 
structural  dynamics  analysis  may  be  assessed  without  the 
necessity  of  dealing  with  aerodynamics  analysis  issues. 

Of  course,  aside  from  the  analysis  itself,  the  results  are 
subject  to  the  accuracy  of  the  experimentally  measured 
airloads  and  blade  loads  as  well  as  the  accuracy  of  the 
physical  properties  of  the  blade  structure.  Other  impor¬ 
tant  considerations  for  this  problem  include  choices  for 
modeling  the  actual  physical  structure,  the  choice  of  de¬ 
grees  of  freedom  are  to  be  included  or  excluded,  and  how 
to  properly  represent  the  boundary  conditions  to  define 
a  well  posed  problem.  In  addition,  some  consideration 
must  also  be  given  to  the  numerical  solution  procedures 
to  insure  convergence  and  practical  computations.  The 
generality  of  the  RCAS  computer  program  provides  op¬ 
portunities  for  developing  complex  structural  models  and 
applying  various  solution  methodologies  tailored  to  the 
problem  at  hand.  The  following  discussion  will  briefly 
survey  the  approach  used  for  this  problem. 

For  this  initial  investigation,  several  assumptions  were 
made  to  bound  the  scope  of  the  problem.  First,  the  prob¬ 
lem  was  limited  to  a  single  isolated  rotor  blade  rotating 
with  constant  angular  velocity  about  a  point  traveling  in  a 
straight  line  with  constant  velocity.  That  is,  any  dynamic 
effects  of  fuselage  or  engine  drive  train  coupling  with  the 
rotor  are  ignored  including  any  coupling  between  blades 
through  swashplate  flexibility. 

The  RCAS  structural  model  is  schematically  depicted 
in  Fig.  24  and  is  composed  of  a  series  of  elements  includ¬ 
ing  rigid  bars,  linear  rotational  and  translational  springs 
and  dampers,  a  rigid  body  mass,  hinges,  a  slide,  and 
twelve  nonlinear  beam  finite  elements.  The  outboard  tip 
of  the  blade  is  swept,  the  structural  and  aerodynamic  axes 
of  the  blade  are  twisted  non-uniformly  along  the  radius,  a 
series  of  rigid  bars  and  spring  elements  are  used  to  repre¬ 
sent  the  pitch  control  linkage  and  a  slide  element  is  intro¬ 
duced  to  control  blade  collective  and  cyclic  pitch  inputs. 
A  linear  lag  damper  element  is  provided,  however,  for 
the  results  presented,  the  actual  measured  damper  force 
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Figure  24:  UH-60  rotor  blade  structural  model  schematic. 


is  used  as  a  blade  excitation  because  the  actual  physical 
damper  characteristics  are  highly  nonlinear.  The  blade 
root  hinge  is  offset  from  the  center  of  rotation  and  the 
actual  elastomeric  bearing  which  accommodates  simul¬ 
taneous  rotations  in  all  three  axes  is  represented  by  three 
discrete  coincident  hinge  elements  for  flap,  lag,  and  pitch 
rotation  of  the  blade.  Physical  properties  for  all  of  the  el¬ 
ements  of  the  blade  model  were  included  as  appropriate. 

To  introduce  the  aerodynamic  normal  force,  chord 
force,  and  pitching  moment,  the  measured  blade  surface 
pressures  were  integrated  at  nine  radial  locations  to  ob¬ 
tain  the  forces  and  moments  per  unit  length,  and  these 
were  interpolated  and  applied  as  discrete  forces  and  mo¬ 
ments,  referred  to  as  mechanical  applied  loads,  at  27  lo¬ 
cations  distributed  along  the  blade.  The  loads  were  ap¬ 
plied  in  the  appropriate  airfoil  chord  coordinates  and  fol¬ 
lowed  the  blade  orientation  as  it  underwent  rigid  flap, 
lead-lag,  and  pitch  motions  as  well  as  elastic  bending  and 
torsion  deformations.  Similarly,  the  blade  loads,  flap- 
wise  bending,  chordwise  bending  and  torsion  moments 
were  calculated  in  the  local  deformed  blade  coordinate 
systems. 

For  the  solution  process,  the  periodic  mechanical 
applied  loads  representing  aerodynamic  and  lead-lag 
damper  were  applied  to  the  rotating  blade  structure  and 
the  solution  proceeded  either  in  1)  the  time  domain  using 
a  time  integration  procedure  until  a  converged  periodic 


solution  for  the  system  degrees  of  freedom  was  obtained 
or,  2)  a  frequency  domain  solution  procedure  was  applied 
based  on  an  assumed  harmonic  solution.  The  later  pro¬ 
cedure  was  most  commonly  applied  as  the  number  of 
harmonics  of  the  solution  could  be  conveniently  speci¬ 
fied  and  this  automatically  filtered  high  frequency  noise 
from  the  experimental  airloads  data.  Ten  harmonics  were 
retained  for  the  results  presented.  To  facilitate  numeri¬ 
cal  convergence  of  the  solutions,  structural  damping  was 
used  to  damp  initial  transients  and  then  withdrawn  as  a 
function  of  time  as  the  solution  converged.  In  addition 
to  the  small  discrete  dampers  used  to  represent  the  blade 
root  hinge  elastomeric  bearing,  0.02%  structural  damping 
was  introduced  into  the  elastic  blade  structure. 

Representing  aerodynamic  forces  as  fixed  mechanical 
loads,  rather  than  being  motion  dependent  as  in  typical 
coupled  aeroelastic  blade  loads  analyses,  eliminates  the 
inherent  aerodynamic  damping  that  typically  aids  the  nu¬ 
merical  solution  convergence  of  aeroelastic  analyses.  For 
the  present  purely  structural  dynamic  response  analysis, 
the  absence  of  damping  can  present  numerical  conver¬ 
gence  problems,  thus  necessitating  the  approach  adopted 
herein.  For  cases  with  very  small  dampers  or  structural 
damping,  problems  may  arise  if  system  natural  frequen¬ 
cies  fall  close  to  the  N/rev  excitation  frequencies.  For  ex¬ 
ample,  the  rigid  flap  natural  frequency  of  an  articulated 
blade  with  zero  hinge  offset  is  1 .0/rev  and,  if  there  is  no 
damper,  response  will  be  infinite  if  the  flapping  moment 
excitation  is  not  zero.  In  such  a  case,  the  1/rev  aerody¬ 
namic  flap  moments  must  be  identically  zero  and  this 
is  consistent  with  the  observed  finite  flapping  response 
of  actual  blades  with  zero  hinge  offset.  However,  if  ex¬ 
perimentally  measured  airloads  contain  any  small  errors, 
the  calculated  dynamic  response  would  be  infinite.  This 
simply  reflects  the  fact  that  calculating  dynamic  response 
from  experimentally  measured  airloads  will  be  impossi¬ 
ble  or  impractical  if  system  natural  frequencies  of  modes 
with  very  small  damping  are  too  close  to  N/rev  excita¬ 
tions.  In  the  present  case,  the  rigid  flap  mode  frequency  is 
approximately  1 .035/rev  and  this  is  sufficiently  far  from 
1/rev  to  enable  reasonably  accurate  flap  motion  response 
to  be  calculated,  even  though  the  1/rev  aerodynamic  flap 
moment  is  very  small. 

Typical  results  for  blade  loads  are  given  in  Fig.  24  - 
Fig.  28  for  the  blade  flapwise  and  edgewise  bending  mo¬ 
ments  and  torsion  moment  at  several  radial  stations  as  a 
function  of  blade  azimuth. 

The  flight  condition  for  these  results  is  level  flight  at 
0.368  advance  ratio.  Calculated  flapwise  bending  mo¬ 
ments  at  50%  blade  radius  are  compared  with  experimen¬ 
tal  measurements  in  Fig.  25  and  the  results  are  shown  to 
be  reasonable.  Chordwise  bending  moments  are  com- 
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Figure  25:  UH-60  blade  flapwise  bending  moment, 
30%  R.  Calculated  blade  loads,  measured  blade  loads. 
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Figure  26:  UH-60  blade  edgewise  bending  moment, 
11.3%  R. 
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Figure  27:  UH-60  blade  edgewise  bending  moment, 
50%  R. 
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Figure  28:  UF1-60  blade  torsion  moment,  30%  R. 


pared  for  two  radial  stations  in  Figs.  26  and  27.  For  the 
1 1.3%  inboard  station,  the  results  are  in  good  agreement, 
reflecting  the  effective  constraint  of  the  lead-lag  damper 
force  applied  at  the  adjacent  damper  attachment  location 
at  8.1%  blade  radius.  At  the  50%  radial  station  the  agree¬ 
ment  is  only  fair,  as  the  measured  peak  bending  moment 
response  is  not  captured  by  the  calculated  result.  It  is  hy¬ 
pothesized  that  this  may  be  related  to  the  constant  rotor 
speed  assumption  in  the  analysis.  A  comparison  of  tor¬ 
sion  moments  at  the  30%  blade  radius  is  given  in  Fig.  28 
and  the  results  are  in  good  agreement  as  far  as  the  oscilla¬ 
tory  behavior  is  concerned.  There  is  a  significant  differ¬ 
ence  in  the  mean  values  of  the  torsion  moment  that  may 
be  due  to  experimental  measurement  error. 

The  present  results  have  been  calculated  assuming  the 
boundary  conditions  for  the  blade  flap  and  lag  motion  are 
free  motion  at  the  blade  root  hinge.  The  blade  root  pitch 
motion  (collective  and  cyclic)  was  input  at  the  controlled 
slide  element  at  the  base  of  the  pushrod.  Strictly  speak¬ 
ing,  this  pitch  control  input  must  be  adjusted  on  subse¬ 
quent  iterations  to  account  for  the  pushrod  elastic  deflec¬ 
tions  in  order  to  match  the  measured  blade  pitch  input. 
An  alternative  approach  would  be  to  use  the  measured 
blade  pitch  and  lag  hinge  motion  as  prescribed  bound¬ 
ary  conditions  for  calculating  the  blade  load  response.  To 
conclude,  the  initial  results  presented  here  are  encourag¬ 
ing  and  provide  a  measure  of  confidence  in  the  capabil¬ 
ities  of  an  advanced  structural  dynamics  computer  pro¬ 
gram  when  applied  to  a  typically  complex  real-world  ro¬ 
tor  blade  configuration. 

Conclusions 

The  hybrid  element  addressing  both  rigid  and  flexible 
body  kinematics  has  been  shown  to  be  effective.  The 
RCAS  implementation  of  the  element  is  consistent  with 
expectations  for  the  cases  investigated.  Error  predic¬ 
tions  for  the  element  correlate  well  with  results.  An  as¬ 
semblage  of  elements  each  with  moderate  deformations 
can  accommodate  large  rigid  body  motions  and  deforma¬ 
tions.  The  tests  exercising  static  energy,  fixed  frame  ki¬ 
netic  energy,  rotating  frame  kinetic  energy,  and  periodic 
boundary  work  produce  consistent  results. 

Furthermore,  the  combination  of  improved  finite  ele¬ 
ment  methods  and  novel  approaches  for  eliminating  sys¬ 
tematic  discrepancies  has  resulted  in  significant  correla¬ 
tion  improvements  with  benchmark  analytical  and  exper¬ 
imental  data.  The  RCAS  nonlinear  beam  element  formu¬ 
lation  is  shown  to  be  adequate  for  beam  elastic  deforma¬ 
tions  beyond  the  range  necessary  for  rotorcraft  applica¬ 
tions.  The  code  is  shown  to  be  applicable  for  complex 


rotor  blade  configurations  as  well. 
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